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ABSTRACT 

Hclioseismic observations have detected small temporal variations of the rotation rate below the solar 
surface corresponding to the so-called 'torsional oscillations' known from Doppler measurements of the 
surface. These appear as bands of slower and faster than average rotation moving equatorward. Here 
we establish, using complementary helioseismic observations over four years from the GONG network 
and from the MDI instrument on board SOHO, that the banded flows are not merely a near-surface 
phenomenon: rather they extend downward at least 60 Mm (some 8% of the total solar radius) and thus 
are evident over a significant fraction of the nearly 200 Mm depth of the solar convection zone. 

Subject headings: Sun: interior — Sun: oscillations — Sun: rotation 



1. INTRODUCTION 

The intensely turbulent state of the solar convection 
zone is revealed by the patterns of granulation, mesogran- 
ulation and supergranulation evident in its surface layers 
(e.g., Brummell, Cattaneo & Toomre 1995). Yet accom- 
panying such turbulent and seemingly chaotic small-scale 
dynamics are also signs of ordered large-scale behavior. 
Most notably the solar differential rotation involves a rel- 
atively smooth decrease in angular velocity from equator 
to pole, both in the surface layers (e.g., Snodgrass 1984) 
and within the convection zone as inferred from helioseis- 
mic measurements (e.g., Thompson et al. 1996; Schou et 
al. 1998a). On the largest scales, the magnetic activity 
similarly exhibits well-defined rules as the 22-year cycle 
progresses. 

An enticing link between the latitudes of field emergence 
and small variations in the rotation rate of the surface lay- 
ers is provided by bands of slightly faster and slower than 
average zonal flows, called torsional oscillations, that were 
observed from direct Doppler measurements to migrate 
towards the equator in a manner similar to the zones of 
solar activity (e.g., Howard & LaBonte 1980; Snodgrass, 
Howard & Webster 1985; Ulrich 1998). Helioseismic anal- 
ysis of data from the Michelson Doppler Imager (MDI) 
instrument (e.g., Scherrer et al. 1995) on the Solar and He- 
liosphcric Observatory (SOHO) spacecraft has confirmed 
the presence of such bands of weak zonal flow, and their 
drift towards the equator, for the present solar cycle (Koso- 
vichev & Schou 1997; Schou et al. 1998a,b; Schou 1999). 
Although the causal relation between these banded flows 
and the zones of magnetic activity is still unclear, it is im- 
portant to understand whether the flows are confined to 
the layer of rotational shear just below the solar surface. In 
this letter, we address such questions using two extensive 
helioseismic data sets, covering slightly over four years, 
obtained with MDI and with the ground-based Global Os- 



cillation Network Group (GONG) project (e.g., Harvey et 
al. 1996). We establish the consistency of the independent 
determinations of the flow from the two data sets, and in- 
fer that the zonal banding signature extends to depths of 
about 60 Mm (or about 8% in radius) below the solar sur- 
face. Thus these are not superficial features, and provide 
evidence of ordered rotational responses as the magnetic 
cycle is progressing. More extensive accounts of such anal- 
yses of zonal flows are provided for GONG data by Howe, 
Komm & Hill (2000) and for MDI data by Toomre et al. 
(2000). 



2. OBSERVATIONS AND DATA ANALYSES 

The rotation rate of the solar convection zone has been 
inferred through inversion of observed rotational splittings 
of solar f and p modes. Two sets of observations have been 
used. One was obtained by the GONG network over the 
period 1995 May 7 to 1999 June 26. This set consists of 40 
overlapping series of 108 days each, with starting dates 36 
days apart. The second set was obtained by MDI and con- 
sists of 11 contiguous sets from 1996 May 1 to 1998 June 
24 (before control of the spacecraft was lost), one set from 
1998 October 23 to 1998 December 21 (after control was 
reasserted), and four sets from 1999 February 3 to 1999 
November 17. While each of the sets nominally covers 72 
days, the ones surrounding the loss of contact are missing 
a few days. 

The dependence of the frequencies on the azimuthal or- 
der was represented in terms of an expansion on orthogo- 
nal polynomials (Ritzwoller & Lavely 1991), expressed in 
terms of the so-called a coefficients a,k(n, I), depending on 
the radial order n and the degree / of the mode, as well as 
on the order k of the coefficient. The odd coefficients are 
related to the angular velocity 0(r, 9) (as a function of the 
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Fig. 1. — Variation of rotation rate with latitude and time from which a temporal average has been subtracted to reveal the migrating 
banded zonal flows, based on OLA inversion of MDI data for target radii of 0.99-R, 0.92/?, 0.88i? and 0.84i?. The rendition is smoothed in 
the temporal direction over a window of 72 days, equal to the period of each separate set of observations. Uniform olive-green vertical bands 
indicate that no data were available in that time period for our analysis. The color bar indicates the dynamic range in nHz of the angular 
velocity. 



distance r to the solar center and the co-latitude 9) by 



27r a 2s +i(n, I) 



K { n f s (r,9) Q(r,9) rdrd9 , (1) 



where R is the solar radius and the kernels K^? are as- 
sumed known from a solar model. The GONG data com- 
prised around 10,000 coefficients, up to a±5, for a total 
of typically 1,200 p-mode multiplets (n,l) for I < 150, 
whereas the MDI sets contained approximately 30,000 co- 
efficients, up to a35, f° r roughly 1,800 multiplets with 
I < 300. The inversions of the relations (1) were car- 
ried out by means of two methods, described by Schou 
et al. (1998a): two-dimensional regularized least-squares 
fitting (RLS) and two-dimensional subtractive optimally 
localized averages (OLA). 

3. RESULTS AND DISCUSSION 

The overall features of the inferred rotation profile are 
very similar to those obtained by Schou et al. (1998a) 
based on 144 days of MDI data. Here we concentrate on 
the time-dependent aspects of the dynamics of the upper 
portions of the convection zone. These are most readily 
studied by considering departures of the ^reconstructed ro- 
tation rate from its temporal average Cl(r, 6). Figure [j] 
shows the evolution of these residuals as a function of lat- 
itude at four target depths, using OLA inversion of the 
MDI data. The inversion is only sensitive to the com- 
ponent of rotation symmetric around the equator; even 
so, to show more clearly the evolution of the features, we 
have included both hemispheres in the plots. The residuals 
show alternating bands of positive and negative zonal ve- 
locity, relative to the temporal mean, converging towards 
the equator with increasing time. The amplitude in these 
flows is around 1.5 nHz, corresponding to velocities of up 
to around 6ms _1 . The flows are visible, at roughly the 
same amplitude, in the inversion targeted at 0.92i?, and 
faint traces are visible in the inversion targeted at 0.88i?; 
we return to the significance of this below. 

Although the signal shown in Figure |l| seems strong and 
coherent, some doubt about its reality may remain. Thus 
access to the independent GONG dataset is essential, with 
the additional advantage that it starts almost a year before 
the MDI data. As a further test, we also apply two differ- 
ent analysis techniques to the set, as shown in Figure |[ 



The GONG and MDI data, where they overlap in time, 
are essentially consistent at r = 0.99i?. At r = 0.95R the 
GONG RLS reconstructions are somewhat noisier so that 
the subtle signature of the migrating zonal bands is less 
obvious. The RLS and OLA inversion results for the MDI 
data agree very well at both depths illustrated. 

To provide a more quantitative comparison, Figure |^ 
shows the GONG and MDI solutions at selected radii and 
latitudes as a function of time. It is evident that the three 
sets of results largely agree within their error bars. Also, 
the variations are highly statistically significant: at the 
equator the angular velocity decreases uniformly, whereas 
at latitude 30°, for example, the variation reflects mainly 
the passage of the band of more rapid rotation. The in- 
versions reveal a considerably different behavior at the two 
depths at latitude 60° , much as could also be inferred from 
Figure 0. 

The interpretation of any inversion results must take 
into account their finite resolution, as well as the proper- 
ties of errors of the inferences. In particular, the solution 
obtained at a given location contains contributions from Q 
at other points, while error correlation between the solu- 
tions at different locations may give the impression of co- 
herent structure where none exists. Quantitative measures 
of these effects can be obtained from detailed analyses of 
the inversion (e.g. Schou et al. 1998a; Howe & Thomp- 
son 1996). As an alternative, we have considered artificial 
data for a number of prescribed rotation laws, with error 
properties corresponding to those of the solar data. The 
mode set and errors used correspond to those of a typical 
MDI set. Superimposed on a smoothly varying flow in the 
latitudinal direction, the artificial rotation profiles possess 
a single pair of flows, 10° wide, and rotating 3 nHz faster 
than the background level, moving equatorwards at 5° per 
sample. There are nine samples in all in each 'time' se- 
quence. Three cases are illustrated in Figure ^ using OLA 
inversions: in the first, the flows extend from the surface to 
a depth of 5% of the solar radius, in the second to a depth 
of 8% and in the third to a depth of 20%. For all cases, 
the branches of the flow are strongly visible at 0.99-R, and 
the flows are evident to about the depth to which they are 
imposed, while disappearing below that depth. The 0.92i? 
case most resembles the solar observations illustrated in 
Figure |l|, whereas the disappearance of the flows in the 
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0.95i? case at 0.93J? and the visibility of the 0.80i? flow 
at 0.847? are both inconsistent with the solar observations. 
This evidence, taken together with many other cases which 
we have tested, strongly suggests that the solar flow struc- 
ture extends to a depth of at least 0.08i? with substantial 
amplitude, but does not extend much further than O.lOi?. 
The latitudinal width of the Sun's banded flows (Fig. [I]) 
is similar to the 10° width assumed in the artificial data: 
this is consistent also with the Doppler measurements at 
the surface (Ulrich, 1998). Comparison of Figures and 
U suggests that the solar flows may be somewhat weaker 
than the 3 nHz chosen for the artificial case. 
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Fig. 2. — Comparison of various inversions of GONG, and MDI 
data showing evolution of residual rotation rate (cf. Fig. hi) at target 
radii 0.99R (left column) and 0.95-R (right column). From the top, 
the three rows show reconstructions with RLS using GONG data, 
with RLS using MDI data, and with OLA using MDI data. 

At the highest latitudes, Figures [l] and || show more dy- 
namical variations than at lower latitudes. This is likely 
related to the lesser moment of inertia associated with the 
polar regions. There may be evidence for the formation 
of a new band of rapid rotation tending towards lower 
latitudes; the further evolution of this feature should be 
followed over the coming years. 
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Fig. 3. — Evolution with time in the residual rotation rate at 
target radii 0.99i? and 0.95R and latitudes 0°, 30°, and 60°. Dia- 
monds (black) show results of RLS inversion applied to GONG data, 
whereas circles (green) and triangles (orange) are for RLS and OLA 
inversions on MDI data. The error bars show one standard devia- 
tion of the inference, as determined from the errors in the observed 
data. 



4. CONCLUSIONS 

Analysis of extended series of GONG and MDI data 
has revealed coherent banded flow structures in the solar 
convection zone. These correspond to the torsional oscil- 
lations detected in direct Doppler observations of the solar 
surface. We have demonstrated that the flows are likely 
to extend to a depth of at least 60 Mm, a substantial frac- 
tion of the total 200 Mm depth of the convection zone, 
and considerably more than the depth (about 35 Mm) at 
which the rotation rate attains its maximum in the sub- 
surface radial shear layer at low latitudes (cf. Schou et 
al. 1998a). In addition, there appear to be other system- 
atic variations with time of the residual rotation rate, with 
different signatures at low and high latitudes (cf. Fig. ||). 

Inversions of global oscillation frequency splittings sam- 
ple the component of rotation symmetric around the equa- 
tor. The actual flows will exhibit some level of asymmetry 
which will depend on the time scale used in the analy- 
sis. Indeed, local analyses by means of the time-distance 
and ring-diagram techniques (Giles, Duvall, & Scherrer 
1998; Haber et al. 2000) have shown features similar to 
those found here, but with clear differences between the 
two hemispheres. 
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Fig. 4. — Sample artificial data tests to study residual rotation 
rate inferred from OLA inversions at various selected radii. The im- 
posed flow used in the calculation of the frequency splittings extends 
from the surface to 0.95-R in the top row, to 0.92i? in the middle 
row, and to 0.80-R in the bottom row. 

The link between the evolving latitudinal positioning of 
the faster zonal bands and of the sites of sunspot emer- 
gence suggest that the dynamics are related, yet how this 
is accomplished is uncertain. The strong magnetic fields 
most likely originate from deep within the Sun, probably 
formed by dynamo action near the base of the convec- 
tion zone (e.g., Spiegel & Zahn 1992; Parker 1993; Weiss 
1994; Charbonneau & MacGregor 1997). Field bundles 
ascending from this region through the convection zone, 
before erupting into the atmosphere as large-scale mag- 
netic loops, could well lead to significant perturbations in 
velocity and thermal fields there. This is likely to be ac- 
companied by some redistribution of angular momentum, 
given that the magnetic structures will attempt to con- 
serve their original angular momentum (e.g., Brummell, 
Cattaneo & Toomre 1995). The coupling of a highly tur- 



bulent medium with ascending magnetic structures, and 
their mutual feedbacks, have not yet been assessed in re- 
cent flux-tube models. Global simulations of turbulent 
convection in rotating spherical shells (e.g., Elliott et al. 
2000; Miesch et al. 2000; Miesch 2000) to study the result- 
ing differential rotation have revealed intrinsic variability 
in zonal flows over intervals of several rotation periods, 
some of which may be inertial oscillations (e.g., Gunther & 
Gilman 1985), but such modelling has not included large- 
scale magnetic fields. Obtaining propagating bands and 
time scales of variation of order the solar cycle seems prob- 
lematic unless there is some selective coupling to magnetic 
processes. Adding to the puzzle is that the evolving zonal 
bands are present at the higher latitudes even before the 
prominent large-scale magnetic eruptions begin (e.g. Ul- 
rich 1998), as within this cycle. Continued helioseismic ob- 
servations as this magnetic cycle is proceeding may help to 
provide clues about such aspects of solar internal dynam- 
ics, for we now have the ability to probe hitherto unseen 
flows well below the solar surface. 
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